SCIENTIFIC 

REPORTS 




OPEN 



SUBJECT AREAS: 

MECHANISMS OF 
DISEASE 

PREDICTIVE MARKERS 



Received 
3 February 2014 

Accepted 
20 March 2014 

Published 
9 April 2014 



Correspondence and 
requests for materials 
should be addressed to 
D.K. (dimitrios- 
karamichos@ouhsc. 

edu] 



In vitro model suggests oxidative stress 
involved in keratoconus disease 

D. Karamichos 1 , A. E. K. Hutcheon 2 , C. B. Rich 4 , V. Trinkaus-Randall 4 , J. M. Asara 3 & J. D. Zieske 2 

1 Department of Ophthalmology, University of Oklahoma Health Sciences Center, USA, 2 Schepens Eye Research Institute/ 
Massachusetts Eye and Ear and the Department of Ophthalmology Harvard Medical School, Boston, MA, USA, 3 Division of Signal 
Transduction/Mass Spectrometry Core, Beth Israel Deaconess Medical Center, Boston, Massachusetts, USA, departments of 
Biochemistry and Ophthalmology, Boston University School of Medicine, 80 E Concord Street, Boston, MA 021 1 8, USA. 

Keratoconus (KC) affects 1 : 2000 people and is a disorder where cornea thins and assumes a conical shape. 
Advanced KC requires surgery to maintain vision. The role of oxidative stress in KC remains unclear. We 
aimed to identify oxidative stress levels between human corneal keratocytes (HCKs), fibroblasts (HCFs) and 
keratoconus cells (HKCs). Cells were cultured in 2D and 3D systems. Vitamin C (VitC) and TGF-03 (T3) 
were used for 4 weeks to stimulate self-assembled extracellular matrix (ECM). No T3 used as controls. 
Samples were analyzed using qRT-PCR and metabolomics. qRT-PCR data showed low levels of collagen I 
and V, as well as keratocan for HKCs, indicating differentiation to a myofibroblast phenotype. Collagen type 
III, a marker for fibrosis, was up regulated in HKCs. We robustly detected more than 150 metabolites of the 
targeted 250 by LC-MS/MS per condition and among those metabolites several were related to oxidative 
stress. Lactate levels, lactate/malate and lactate/pyruvate ratios were elevated in HKCs, while arginine and 
glutathione/oxidized glutathione ratio were reduced. Similar patterns found in both 2D and 3D. Our data 
shows that fibroblasts exhibit enhanced oxidative stress compared to keratocytes. Furthermore the HKC 
cells exhibit the greatest level suggesting they may have a myofibroblast phenotype. 



Keratoconus (KC) is a progressive degenerative disease and is a major clinical problem worldwide. KC affects 
4 to 600 per 100,000 people and in advanced cases can lead to severe loss of vision. KC affects the cornea 
structurally, causing it to thin, weaken and protrude into a more conical shape. While KC appears in the 
teenage years to young adulthood, its rate of progression is generally unpredictable. Currently, the use of contact 
lenses is used at the early stages of the disease, whereas advanced stages often require surgery, such as corneal 
transplantation, to maintain or improve vision. Surgeries, however, do not come without complications, which 
include vascularization of the corneal tissue and rejection of the donor cornea. Corneal collagen crosslinking 
(CXL) is a new upcoming treatment for keratoconus disease where riboflavin is added to the cornea and is 
activated by ultraviolet (UV) light 1,2 . Riboflavin has been shown to increase collagen crosslinking in the cornea 
and increase its strength leading to less deformation. The pathogenesis of KC is still unclear, however, a variety of 
factors have been considered, including genetics and cellular mechanisms, as well as oxidative stress 3 " 9 . One of the 
few recognized links to KC disease is the association of individuals with Down syndrome. Approximately, 15% of 
the individuals with Down syndrome exhibit KC disease 6,10 . The common gene defect between the two is the 
superoxide dismutase-1 (SOD1) where it has been proposed as a possible candidate gene for familial KC 11 . Other 
gene defects associated with KC disease include VSX1, COL6A1, COL8A1, and MMP9 as reviewed by Nowakand 
Gajecka 12 . 

Oxidative stress is an indication of the accumulation of reactive oxygen species (ROS) that may lead to 
disturbances in the normal redox state of the cell and cause a toxic effect by free radicals and damage all cellular 
components, including DNA, proteins and lipids. In humans, oxidative stress has been linked to a variety of 
diseases, such as cancer 13 , Parkinson's disease 14 , heart failure 15 , and myocardial infraction 16,17 . In KC, the mRNA 
levels of antioxidant enzymes have been investigated 18 in KC cornea buttons. The authors concluded that KC 
corneas undergo oxidative stress and tissue degradation. In agreement with this study, are several others where 
the antioxidant activities as well as protein levels have been investigated 19,20 . 

Previously, we used a novel 3D model to study KC, which utilizes human keratoconus cells (HKCs) incubated 
in the presence of a stabilized derivative of Vitamin C (VitC) and TGF-fj3 (T3). Under these conditions the HKCs 
deposited an extracellular matrix (ECM) that resembles normal stromal ECM 21 . In that study HKCs without T3 
secreted an ECM with fibrotic characteristics. 
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Figure 1 | Representative image showing the morphology of HKCs in 3D 
cultures. HKCs maintained fibroblastic, elongated morphology under all 
conditions. Red = Phalloidin, Blue = TOPRO-3. Bars = 50 microns. 

Targeted metabolomics using mass spectrometry has been used to 
profile many diseases such as cancers from biological tissue sources 
such as cerebral spinal fluid (CSF), plasma and tumor tissue 22 " 25 . In 
our current study, we investigated the metabolic differences between 
human corneal keratocytes (HCKs), fibroblasts (HCFs) and HKCs, 
under conventional 2D cultures and using our 3D in vitro model. 
Investigating both systems allow us to compare cells in a monolayer 
culture and cells in a self-assembled ECM, which more resembles an 
in vivo like condition. From these studies, our data indicated that 
HKCs express metabolites that are indicative of oxidative stress both 
in the 2D and 3D cultures. 

Results 

Cell morphology and real time-PCR. We investigated the morpho- 
logy of the three cell types. Both HKCs and HCFs were elongated or 
fibroblastic in appearance, while the HCKs were dendritic in shape. 
Figure 1 shows the morphological phenotype of HKCs which is very 
similar to HCFs as we have shown previously 28 . In 3-D cultures, 
HCKs maintained their dendritic morphology and were distinc- 
tively different from HCFs as previously shown 28 . The morpholo- 
gical characteristics were maintained both with cells in culture and 
cells on polycarbonate membranes. We did not notice any significant 
morphological changes between initial cell expansion and experi- 
mental conditions. 

We further characterized the cells by investigating several essential 
probes with qRT-PCR for corneal stromal cells: Col I, III and V. Col I 
and V are routinely found in healthy adult corneal stroma, whereas 
Col III is a sign of myofibroblast differentiation and a wounded/ 
scarred cornea. We have shown previously 28 that Col I and V (~4- 
5 fold) are significantly upregulated in HCKs and HCFs when T3 is 
present, independent of VitC stimulation 28 . In Table 1 we compared 
the two cell types to HKCs. Expression levels of Col I and Col V, for 
HKCs, were identical to HCFs; however a 4 fold decrease was seen 
when compared to HCKs for both Col I and -V. Col III (Table 1) 
expression was significantly higher (~4 fold) with HKCs when com- 
pared to HCFs and even more so when compared to HCKs ( — 108 
fold). Since Col III is a marker for corneal fibrosis and myofibroblast 
population presence, this is indicative of the different cellular state of 
the KC corneas. Our data suggest that HKC are myofibroblasts in 
agreement to previous study where high smooth muscle actin (SMA) 
expression was shown 21 . 



Table 1 qRT-PCR analysis for all three cell types HCKs, HCFs, and 
HKCs for Collagen-I, -III and -V. Col I and Col V expression was 
down-regulated 4 fold in HCFs and HKCs when compared to HCKs. 
Col III was significantly higher (~4 fold) with HKCs when compared 
to HCFs and 1 08 fold compared to HCKs 

Collagen-I Collagen-Ill Collagen-V 

(RQ Mean) (RQ Mean) (RQ Mean) 

HCK 11.46 + 0.68 0.07667 ± 0.001 6 1 1 .43 ± 0.2706 
HCF 2.84 ± 0.30 2.27 ± 0.15 1.773 ±0.1877 

HKC 2.763 ± 0.29 8.637 ± 0.61 3.12 ±0.26 



Predicted Pathways of Regulation. Our system allows us to investi- 
gate the levels of various metabolites, both in cells that are conven- 
tionally cultured (2D), as well as in self-assembled 3D ECM. From 
our overall raw data analysis, we observed that out of 249 targeted 
metabolites, the number of those robustly detected ranged from 127 
to 155 across the samples with a mean value of 144. In Figure 2 and 3 
the predicted metabolic pathways that are affected based on the raw 
data is shown in 2D (Figure 2) and 3D (Figure 3) systems. When we 
analyzed the metabolites that were up regulated by at least 2 -fold in 
HKCs compared to HCKs a number of pathways that are involved in 
oxidative stress were affected. Examples are the citric acid cycle, the 
glycine-serine-threonine, the betaine and the glutathione metabo- 
lism. On the other hand, based on the metabolites where HCKs 
expressed in higher levels than HKCs, our software indicated 
significant regulation of pathways that are mainly involved in cell 
growth (biotin metabolism), toxicity protection (ammonia metabo- 
lism, and apoptosis prevention (ubiquinone biosynthesis). Based on 
these results we further analyzed the metabolites involved in 
oxidative stress. 

Lactate: Cells and Constructs. Our system allows us to investigate 
the levels of various metabolites, both in cells that are conventionally 
cultured (2D), as well as in self-assembled 3D ECM. Lactate is 
produced from pyruvate via glycolysis and is an important 
metabolite that has been implicated in corneal edema 31 , When 
lactate levels are high, corneal edema appears. High lactate levels 
are also a hallmark of many diseases and lactate has become a 
potential target for cancer therapies 32 . Figure 4 shows a clear 
difference in lactate levels between cells cultured in monolayer 
(2D) and stratified cultures (3D). Lactate levels were significantly 
(p<0.05) higher when cells were allowed to develop their own 
ECM and grow in an environment that mimics the in vivo. HKC 
were significantly higher than both HCFs and HCKs. T3 stimulation 
only significantly affected HCFs (p<0.05) with a 2 fold up regulation 
of lactate levels. HKCs showed down regulation of lactate (p<0.05) 
with the same T3 treatment. Similarly, in the monolayer system, 
HCKs maintained the lowest levels of lactate when compared to 
both HCFs and HKCs. 

Lactate/Malate ratio: Cells and Constructs. Lactate has been linked 
to oxidative stress and when paired with other specific metabolites, 
can act as an indicator of oxidative stress 33 . One such combination is 
the lactate-malate ratio (L/M), which is known to be a critical 
regulator of oxidative stress 34,35 . Malate is an intermediate in the 
citric acid cycle (or TCA cycle) produced from fumerate and an 
indicator of citric acid cycling. In vivo, the degree of elevation of 
the ratio and its timing are related to the severity of hypoxia and 
oxidative stress 35 . In our systems, HCKs showed the lowest L/M ratio 
independent of condition and treatment (Figure 5), when compared 
to both HCFs (p<0.05) and HKCs (p<0.01). In agreement with our 
lactate data, HKCs seem to be under more oxidative stress. This is 
enhanced in the 3D system where there is a significant 4 fold 
difference from the HCF (p<0.01). When HCFs were treated with 
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Figure 2 | Summary of pathway enrichment analysis in 2D system. Above is a display of the diversity of signaling pathways that are enriched on the basis 
of all the metabolites passing filtering criteria. The most significant p-values are in red while the least significant are in yellow and white. A) Shows 
pathways affected based on metabolites that were at least 2-fold down-regulated in HKCs compared to HCKs, B) Shows pathways affected based on 
metabolites that were at least 2-fold up-regulated in HKCs compared to HCKs. 



T3, L/M ratio was negatively affected showing a 2 and 4 fold up 
regulation in 2 (p<0.05) and 3D system (p<0.01) respectively. 17 
M ratio, on the other hand, was reduced in HKCs (p<0.05). Overall, 
HKCs abnormal levels of oxidative stress are enhanced in the 3D 
cultures, which is a system more accurately mirroring in vivo 
regulations (Figure 5). Lactate/malate can indicate the relative 
levels of glycolytic activity to TCA cycle activity. 

Lactate/Pyruvate ratio: Cells and Constructs. Pyruvate represents 
the final step of glycolysis and a precursor to lactate. The relative ratio 



of lactate/pyruvate (L/P) has been used as a measure of anaerobic and 
aerobic corneal metabolism. However, the ratio also has been 
proposed as an indicator of oxidative stress 28 . An increase in the L/ 
P ratio has been reported to accompany corneal degradation and 
oxidative stress 33 . Figure 6 shows ratio regulation on both our 
systems for all cell types. HCKs showed the lowest ratio of L/P in 
both systems and all conditions (p<0.001). HKC on the other hand 
displayed the greatest up regulation of L/P (p<0.0001) in the 3D 
model indicating oxidative stress. Upon T3 stimulation L/P levels 
were positively affected and significantly down regulated (p<0.001). 
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Figure 3 | Summary of pathway enrichment analysis in 3D system. Above is a display of the diversity of signaling pathways that are enriched on the basis 
of all the metabolites passing filtering criteria. The most significant p-values are in red while the least significant are in yellow and white. A) Shows 
pathways affected based on metabolites that were at least 2-fold down-regulated in HKCs compared to HCKs, B) Shows pathways affected based on 
metabolites that were at least 2-fold up-regulated in HKCs compared to HCKs. 
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Figure 4 | Lactate metabolite regulation in 2D and 3D culture models for 
HCKs, HCFs, and HKCs. Two conditions were examined and analyzed for 
all cell types: Control: VitC only, and T3:VitC + TGF-fi3. Lactate levels 
were significantly (p<0.05) higher in 3D systems compared to 2D. Upon 
T3 stimulation only HCFs showed significant up-regulation of Lactate 
levels (p<0.05). 

Consistent with our lactate results, the L/P levels were enhanced in 
the 3D model for the HCKs. No differences were seen for the HCF 
which showed identical regulation on all systems and conditions. 

GSH ratios: Cells and Constructs. Another way to measure levels of 
oxidative stress is by looking at a highly studied antioxidant and 
metabolite glutathione. Glutathione exists in both reduced (GSH) 
and oxidized (GSSG) states. A decreased GSH-to-GSSG ratio is 
considered indicative of oxidative stress 36,37 . In the keratoconic 
cornea, levels of GSH are low in central corneal buttons 38 . In this 
study, we investigated the level of GSH-to-GSSG ratio in our 2D and 
3D systems. (Figure 7) The lowest ratio of GSH-to-GSSG was 
detected in the HKCs (p<0.05) under all systems and conditions, 
indicating that the cells were and remained under oxidative stress. 
On the other hand the highest ratio, not surprisingly, was seen with 
HCKs (p<0.05) indicating less oxidative stress levels. Unlike other 
indicators so far, such as lactate, L/M and L/P ratios (Figures 4, 5, and 
6 respectivelly), the 3D system seemed to recover some of the HKCs 
stress levels, i.e under less oxidative stress. Figure 7 shows higher 
values of GSH-to-GSSG ratio for all three cell types in the 3D 
system compared to 2D (p<0.05). T3 treatment significantly 
increased the ratio values for HCKs (p<0.05) and HCFs (p<0.05), 
but not for HKCs. Overall, our data throughout showed that HKCs 
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Figure 6 | Lactate to Pyruvate metabolite ratio (L/P) in 2D and 3D culture 
models for HCKs, HCFs, and HKCs. Two conditions were examined and 
analyzed for all cell types: Control: VitC only, and T3:VitC + TGF-|33. 
HCKs showed the lowest ratio of L/P under all conditions (p < 0.00 1 ) . HKC 
displayed the greatest up regulation of L/P (p<0.0001) in the 3D model. T3 
stimulation reversed that (p<0.001). 

are distinct from both HCKs and HCFs and are under more oxidative 
stress regardless of which system is used (2D or 3D). 

Arginine: Cells and Constructs. Arginine is another important 
metabolite that has been linked to oxidative stress. In the cornea, 
arginine is a key player in the overall immune privilege of the eye. 
Inhibition of arginine has been shown to accelerate graft rejection 35 . 
We investigated the regulation of arginine in our systems and the 
results are shown in Figure 8. As expected, HCKs arginine levels were 
significantly higher than HCFs and HKC (p<0.01), however that 
was only for the 2D system. In the 3D model, HCKs and HKCs 
showed similar levels where HCFs were significantly lower. This is 
another metabolite where its expression is regulated differently when 
the cells secrete their own ECM versus the conventional cell culture 
2D model. While there is a clear indication that HKCs' arginine levels 
are lower than the native HCKs, and therefore under oxidative stress, 
our data with T3 is somewhat confusing. In 2D, HCKs showed a 
massive down regulation of arginine, where HCFs and HKCs did not; 
where in 3D the effect was inverted with a big up regulation of 
arginine expression in HCKs. Clearly the role of T3 in oxidative 
stress regulation is not well understood and further characteriza- 
tion is currently under way. 
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Figure 5 | Lactate to Malate metabolite ratio (L/M) in 2D and 3D culture 
models for HCKs, HCFs, and HKCs. Two conditions were examined and 
analyzed for all cell types: ControhVitC only, and T3:VitC + TGF-|33. 
HCKs showed the lowest L/M ratio under all conditions compared to both 
HCFs (p<0.05) and HKCs (p<0.01). T3 significantly increased L/M ratio 
in HCFs (p<0.01) and decreased it in HKCs (p<0.01). 
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Figure 7 | Glutathione reduced (GSH) to Glutathione oxidized (GSSG) 
ratio (GSH-to-GSSG) in 2D and 3D culture models for HCKs, HCFs, and 
HKCs. Two conditions were examined and analyzed for all cell types: 
Control: VitC only, and T3:VitC + TGF-|33. HKCs showed the lowest ratio 
under all conditions (p<0.05). T3 treatment significantly increased the 
ratio values for HCKs (p<0.05) and HCFs (p<0.05), but not for HKCs. 
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Figure 8 | Arginine metabolite regulation in 2D and 3D culture models 
for HCKs, HCFs, and HKCs. Two conditions were examined and analyzed 
for all cell types: Control : VitC only, and T3:VitC + TGF-F33. HCKs 
arginine levels were significantly higher than HCFs and HKC (p<0.01) 
only for the 2D system. 

Discussion 

Despite the huge advancements in technology, very little is known 
about metabolic activity of corneal cells derived from Keratoconus 
patients. This is surprising considering that metabolomics has been 
used successfully in ocular diseases 19,20,3139,40 . We therefore investi- 
gated the metabolic activity of corneal cells derived from both normal 
individual as well as Keratoconus patients. In order to more accur- 
ately quantify metabolic activity and correlate it to what is happening 
in vivo we used two in vitro models; a) conventional 2D cultures and 
b) a 3D in vitro model that we previously established 27 as a novel way 
to study keratoconus disease at the cellular level in vitro. Our data 
clearly shows that HKCs are significantly different in terms of their 
metabolism in vitro, when compared to corneal fibroblasts. Specific 
metabolites that are known to be key players in oxidative stress were 
up regulated with HKCs as compared to HCFs indicating oxidative 
stress status. In turn, both of these cell types were different from 
normal keratocytes. Almost all metabolites tested here indicated that 
HCKs were under less stress when compared to either HCFs or 
HKCs. 

Oxidative stress is defined as the imbalance between the systemic 
manifestation of reactive oxygen species and a biological system's 
ability to detoxify or repair any resulting damage from the reactive 
intermediates 41 . In humans, oxidative stress depends on the size of 
these changes and whether the cell can overcome them and return to 
its original stage. One of the main sources of reactive oxygen is the 
leakage of activated oxygen from mitochondria during the metabolic 
pathway of oxidative phosphorylation 42 . The role of the mitochon- 
dria in energy metabolism and oxidative stress has been well estab- 
lished 43 45 . It appears to be intimately linked to their generation of 
ROS that play a regulatory role in cellular metabolic processes. In 
fact, the majority of ROS are products of mitochondrial respiration. 
ROS production and oxidative stress have been linked with an array 
of pathologies, including type II diabetes, atherosclerosis, ischemia/ 
reperfusion injury, and others 46-48 . 

In ocular surface there have been multiple reports about oxidative 
stress, including keratoconus. Chwa and co-authors 49 reported that 
KC fibroblasts have an inherent, hypersensitive response to oxidative 
stressors that involves mitochondrial dysfunction and mtDNA 
damage. Nguyen and co-workers 40 recently showed that lactate 
transport was a significant component in the corneal endothelial 
pump and cornea edema in vivo could be avoided if high levels of 
lactate were cleared, this was in agreement with previous studies by 
Klyce 31 . In our system, lactate levels were higher for the Keratoconus 
cells in the 3D model when compared to Fibroblasts and in both 2D 
and 3D systems when compared to keratocytes. Fu and coworkers 
showed that another metabolite, arginine, plays a key role in cornea, 



and the overall immune privilege of the eye in vivo 39 . Arginine inhibi- 
tion led to accelerated corneal graft rejection. Our data shows that 
keratocytes had higher level of arginine than both fibroblasts and 
keratocytes. HCK Control and T3 treated seem to be significantly 
regulated between the two systems. Taking these together it is clear 
that Keratoconus cell metabolism is different from normal corneal 
keratocyte and may provide clues to diagnose and treatment of the 
disease. Interestingly, both lactate and arginine have been previously 
linked to corneal edema and thinning 31,39,40,50 which are characteris- 
tics of Keratoconus disease. 

Another metabolite that we identified with abnormal levels in the 
HKCs was glutathione (GSH), an antioxidant peptide that prevents 
cellular damage from oxidative stress 51 . This is in agreement with 
previous studies showing that GSH contribute to corneal degenera- 
tion in keratoconus 52,53 . In order to ensure accuracy of our results we 
plotted ratios of GSH with other metabolites that are known to be a 
sign of oxidative stress such as malate and pyruvate. In all cases 
HKCs showed such levels suggesting that were under more stress 
than keratocytes. 

The cornea is a transparent avascular tissue that absorbs approxi- 
mately 80% of the incident ultraviolet B (UVB) light 54 , making it 
highly sensitive and vulnerable to damage from free radicals and 
ROS 54,55 . In the healthy cornea, there are a number of defensive 
mechanisms that are present to minimize and reduce the risk of 
oxidative damage 53,56 . A variety of antioxidant enzymes are present 
such as catalase, glutathione peroxidase, and glutathione reductase 
ensuring the removal of ROS generated by UV light exposure 57 . Cells 
will also defend themselves with enzymes such as alpha- 1 -microglo- 
bulin, superoxide dismutases, glutathione peroxidases and peroxi- 
doxins 58,59 . Small molecules antioxidants such as Vitamin C, 
topopherol, and glutathione also play a vital role in the cornea s 
defense mechanism. 

Collectively our data shows that cells derived from Keratoconus 
patients are under oxidative stress both in conventional 2D model as 
well as when they secrete their own ECM (ie 3D model). Our pathway 
prediction software supports this conclusion, as many of the altered 
pathways in HKCs were oxidative stress related. While this is pre- 
dictive software it gives us a good indication when alterations 
between cell types are present. Based on these results future studies 
should help identify and characterize specific pathways that poten- 
tially might lead to metabolic therapies for keratoconus. This is a 
novel way to study the metabolism of these cells and identify the key 
players that can help us understand the mechanism behind this 
disease. Ongoing studies aim to identify more key metabolites as well 
as further characterize the role of TGF-b3 in regulating metabolism. 
Form the data shown here T3 seems to have differential effect on 
multiple metabolites making its role hard to identify. 

Methods 

Primary cultures of Human Corneal Keratocytes (HCK), Fibroblasts (HCF) and 
Keratoconus cells (HKC). HCKs and HCFs were isolated from human corneas from 
healthy patients without ocular disease. All samples were obtained from NDRI 
(National Disease Research Interchange; Philadelphia, PA). HKCs were isolated from 
human corneas from patients with Keratoconus defects. These samples were obtained 
from Ula Jurkunas (Massachusetts Eye and Ear Infirmary, Boston, MA, USA) and 
Jesper Hjortdal (Aarhus University Hospital, Aarhus, Denmark). All research 
adhered to the tenets of the Declaration of Helsinki. Tissue was processed, as 
previously described 26 . Briefly, corneal epithelium and endothelium were removed 
from the stroma by scraping with a razor blade. The stromal tissue was then cut into 
~2 X 2 mm pieces and placed into T25 culture flasks. Explants were allowed to 
adhere to the bottom of the wells and then Eagle's Minimum Essential Medium 
(EMEM: ATCC; Manassas, VA) with either 1% (HCKs) or 10% (HCFs and HKCs) 
fetal bovine serum (FBS: Atlantic Biologicals; Miami, FL) was added. Following 1-2 
weeks of cultivation, the cells were passaged into 100 mm cell culture plates and 
allowed to grow to 100% confluence before being used in our two systems (2D and 
3D). 

Conventional 2D cultures. All cell types (HCKs, HCFs and HKCs) were cultured on 
conventional 6-well tissue culture plates and processed for qRT-PCR. Cells were 
seeded at 10 6 cells/well and cultured in EMEM with either 1% (HCKs) or 10% FBS 
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(HCFs and HKCs) + VitC ± T3 (0.1 ng/ml TGF-{33: R&D systems, Minneapolis, 
MN); (Vit C: 0.5 mM 2-O-a-D-glucopyranosyl-L-ascorbic acid: Wako Chemicals 
USA, Inc.; Richmond, VA). Cultures without T3 served as Controls (C). 

Real time PCR (qRT-PCR). Total RNA was extracted and processed, and qRT-PCR 
was performed, as previously described 21 .Total RNA was extracted from the cells 
using Genejet RNA Purification Kit (Genejet RNA Purification Kit: 
ThermoScientific, K0731). Genomic DNA was removed by incubation with RNase- 
free DNase I (New England BioLabs, M0303S) in the presence of RNase inhibitor. The 
RNA was annealed with oligo dt and random hexamer primers, and first strand 
synthesis was carried out with MuLV reverse transcriptase (Life technologies, Grand 
Island, NY). Negative controls were performed without reverse transcriptase. A Vii7A 
(Life technologies, Grand Island, NY) was used to perform the qRT-PCR with ABI 
TaqMan gene expression assays — CollAl: Hs00164004_ml, Col3Al: 
Hs00943809_ml, Col5Al: Hs00609088_ml, and Keratocan: Hs00559942_ml and 
the Eukaryotic 18S rRNA Endogenous Control, 4308329. Results were calculated 
using the AACt method, using 18S rRNA as the endogenous control. 

Assembly of Extracellular Matrix - 3D constructs. As previously described 21 ' 27 ' 28 , 
cells were plated on transwell 6-well plates containing polycarbonate membrane 
inserts with 0.4 um pores (Costar; Charlotte, NC) at a density of 10 6 cells/ml. Cells 
were cultured in EMEM with either 1% (HCKs) or 10% FBS (HCFs and HKCs) + 
VitC ± T3 (0.1 ng/ml). The cultures were allowed to grow for 4 weeks. The optimal 
concentration of T3 was previously determined by comparing a concentration series 
ranging from 0.1-10 ng/ml 27 . Cultures without T3 served as Controls (C). Constructs 
thicknesses were similar with our previous studies 21 ' 27 28 . 

Metabolite extraction. All tissues were collected and processed as previously 
reported 24,29 . Briefly, samples were grounded using a homogenizer in ice-cold 80% 
MeOH. Samples were centrifuged (14,000 g, 10 min, 4°C) and supernatants 
incubated on dry ice. Pellets were further disrupted in ice-cold 80% MeOH and 
combined with previous supernatants. Plasma metabolites were extracted twice in 
80% ice-cold MeOH. Metabolite extracts were vortexed and centrifuged (14,000 g, 
10 min, 4°C). Supernatants were evaporated and stored at — 80°C until further 
analysis. 

Targeted Mass Spectrometry. Samples were re-suspended using 20 uL HPLC grade 
water for mass spectrometry, and 5-7 ul was injected into a hybrid 5500 QTRAP 
triple quadrupole mass spectrometer (AB/SCIEX) coupled to a Prominence UFLC 
HPLC system (Shimadzu, Columbia, MD) 24,29 . Samples were analyzed via selected 
reaction monitoring (SRM) of a total of 256 endogenous water soluble metabolites for 
steady-state analyses. Some metabolites were targeted in both positive and negative 
ion mode, for a total of 289 SRM transitions, using positive/negative ion polarity 
switching. ESI voltage was +4900 V in positive ion mode and —4500 V in negative 
ion mode. The dwell time was 3 ms per SRM transition and the total cycle time was 
1.55 seconds. Approximately 10-14 data points were acquired per detected 
metabolite. Samples were delivered to the mass spectrometer via hydrophilic 
interaction chromatography (HILIC) using a 4.6 mm i.d X 10 cm Amide XBridge 
column (Waters) at 375 jj.l/min. Gradients were run as follows: 1) 0-5 minutes with 
85 to 42% buffer B (HPLC grade acetonitrile),; 2) 5-16 minutes with 42 to 0% buffer B; 
3) 16-24 minutes with 0% buffer B; 4) 24-25 minutes with 0 to 85% buffer B; and 5) 7 
minutes with 85% buffer B. This last step was to re -equilibrate the column. Peak areas 
from the total ion current for each metabolite SRM transition were integrated using 
MultiQuant v2.0 software (AB/SCIEX). 

Pathway Enrichment Analysis. Pathway enrichment (representation) analysis was 
performed using an online freeware program; Metaboanalyst (www.Metaboanalyst. 
ca) intended for the analysis of metabolomics data. 30 Briefly, we performed the 
analysis using only the metabolites that were up or down regulated by 2-fold as 
indicated by our statistical analysis. We chose 2-fold cutoff in order to ensure that we 
included only the vastly abundant metabolites. Only metabolites that were present on 
all biological samples were considered for further analysis. The metabolites that 
passed our criteria were input into the software and the pathway enrichment analysis 
was run. The output of this algorithm marks a metabolic pathway if the input list of 
the metabolites suggests that the pathway associated with a specific pathway are 
present at levels that would not be expected by random chance. 

Statistical Analysis. All experiments were repeated at least three times (three separate 
donors) and data was analyzed for significant variations (p<0.05) using the Student's 
t-test and Dunnett's Multiple Comparison test. 
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